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Ascertaining the core-collapse supernova mechanism is a complex, and yet unsolved, problem 
dependent on the interaction of general relativity, hydrodynamics, neutrino transport, neutrino- 
matter interactions, and nuclear equations of state and reaction kinetics. Ab initio modeling of 
core-collapse supernovae and their nucleosynthetic outcomes requires care in the coupling and 
approximations of the physical components. We have built our multi-physics Chimera code for 
supernova modeling in 1-, 2-, and 3-D, using ray-by -ray neutrino transport, approximate general 
relativity, and detailed neutrino and nuclear physics. We discuss some early results from our 
current series of exploding 2D simulations and our work to perform computationally tractable 
simulations in 3D using the "Yin- Yang" grid. 
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1. Simulating the revival of the supernova shock 

Reviving the stalled shock above the nascent neutron star in core-collapse supernovae (CC- 
SNe) has been a vexing and compelling problem in computational astrophysics for more than 40 
years. The cumulative analysis of CCSN simulations has shown that there are many physical in- 
gredients that impact the shock revival, some enhancing revival and others suppressing it [1, 2, 3]. 
Important microphysical inputs include the various neutrino-matter interactions and the nuclear 
equation of state (EoS), while macroscopic inputs include the treatment of hydrodynamics, general 
relativity, and neutrino transport. 

Shock revival is subject to physical instabilities that make CCSNe inherently an asymmetric 
problem. Neutrino-driven convection not only induces non-spherical motion, but more critically 
boosts the rate of energy deposition in the 'hot mantle' between the proto-neutron star (PNS) and 
the shock. The shock's instability against non-radial perturbations, the Standing Accretion Shock 
Instability (SASI), is dominated in 2D (axisymmetry) by an £ = 1 'sloshing' mode that expands 
the volume inside the shock. The flow onto the PNS is channelled into (intermittent) accretion 
streams. These features make the dynamics of CCSNe intrinsically non- spherical. Here we report 
on some aspects of our current series of 2D simulations, which show shock expansion to many 
thousand km and 'diagnostic' explosion energies [4] close to the canonical 1 B (= 10^^ erg) [5] 
and look ahead to our forthcoming 3D simulation. The 2D simulations begin from the 12, 15, 20, 
and 25 Mq progenitors of Woosley and Heger [6] with 512 adaptive radial zones and 256 angular 
zones covering the entire sphere in our Chimera code [5]. The fluid motion is solved using 
the VHl implementation of the dimensionally split, Lagrangian-remap (PPMLR) scheme, with a 
multipole expansion of gravity that includes a GR monopole term. The neutrino transport is a "ray- 
by-ray-plus" (RbR) scheme using GR-enhanced, multi-group flux-limited diffusion [7] with the 
best available neutrino-matter interactions. We utilize the K = 220 MeV version of the Lattimer- 
Swesty [8] EoS for p > 10^^ g cm"^ and Cooperstein [9] EoS for NSE where p < 10^^ g cm"^ 
and an extended version of the Cooperstein ^"^^y-EoS throughout. Outside NSE conditions we 
use the integrated 14-species XNet [10] a-network. 

2. Dynamics during shock revival 

After a brief epoch (^10-40 ms after bounce) of prompt convection of lepton and entropy un- 
stable gradients left in the wake of the rebound shock, the hot mantle becomes relatively quiescent 
until neutrino heating results in convective instability. By 90 ms after bounce in our 12-Mq model 
(Fig. 1; top) ^10 convective plumes are rising from the heating/cooling interface (gain radius), 
but have yet to distort the spherical shape of the shock. As the plumes grow in size they reduce 
in number and help trigger growth of the SASI sloshing mode. These instabilities help to form a 
configuration common to 2D CCSN simulations [11, 12, 13, 14] with the shock elongated along the 
axis, large plumes rising near the poles, and an accretion stream forming at the shock cusp between 
lobes where the shock radius is typically at a minimum. A different configuration found in some 
2D CCSN simulations can occur if the accretion stream forms at, or near, one pole resulting in a 
single lobed configuration. This configuration does not appear in this series of models, but it did 
occur in the 20 Mq model in our previous series [15] and some models of [16, 4, 14]. 
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Figure 1: Net neutrino heating (red) and cooling (blue) on a grayscale entropy background for 12 Mq model 
at 90 ms (top) and 240 ms (bottom) after bounce. The transition from heating to cooling is known as the 
gain 'radius'. Darker grays indicate lower entropy regions. 



At 240 ms (Fig. 1 ; bottom) the shock has a minimum /^sh ^350 km near the equator with a low 
entropy stream (darker shades). In the inner 200 km, there are 2 major streams reaching the PNS 
near the equator and on the left and a thinner stream on the right. Where these penetrate to the PNS 
there is a localized, intense heating as shown by the deep red color. Because of the RbR transport, 
the strongest heating is only above the 'hot spots' where the streams decelerate at the PNS. We 
would expect to find some lateral spreading of these intense heating regions in simulations with 
full multi-D transport (i.e., without the RbR-approximation) [17]. Along the right (+) pole where 
there is no active accretion, there is less intensive and smoother heating in a polar outflow. 

A useful approach to diagnosing fluid motion is with 'Lagrangian tracer particles'. These 
'tracers' are initialized in equal-mass spaced 'rows' of 40 tracers with an equal volume lateral 
spacing within a row [18]. Fig. 2 shows the trajectories of a row of tracers that start from inside 
the O-rich Si-layer of the progenitor. This row of tracers strikes the shock when /^sh ^500 km and 
the shock is already significantly distorted and beginning to revive. The tracer closest to each pole 
in this row is deflected into a lobe of the expanding shock of the nascent CCSN explosion. Most 
of the remaining tracers are deflected along the inward-curving shock into the equatorial accretion 
stream and plunge toward the PNS. Tracer analysis provides an alternative view of the supernova 
dynamics than the commonly used 'snapshot' images. We end the plotting of these traces at the 
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Figure 2: Trajectories for a rows of tracers initialized at M = 1. 393 Mq in the Si-shell in the 12 Mq model 
colored by starting angle. 

time when most of the tracers reach the PNS to illustrate in the infall stage of their motion, but 
subsequent tracer analyses will help identify how material fills the plumes and becomes ejecta. 

3. Ejecta 

Though these models are not complete enough for detailed nucleosynthesis (the 12 Mq model 
at 850 ms after bounce with an asymptotic explosion energy and negligible accretion seems to 
be nearing that point [5]), by examining the electron fraction, Ye, of the ejecta at 450 ms for all 
models (Fig. 3) we can see that proton-richness of the ejecta appears to increase with progenitor 
mass. For the 12 Mq model, there is more lower- (lighter blues) moderately neutron-rich ejecta 
than higher-Fg, moderately proton-rich ejecta (red-yellow), and Ye < 0.52. For the 25 Mq model 
there is a shift to higher-F^ ejecta with notable fractions above Ye = 0.52, though there remains 
some low-Ye material visible. The 15 and 20 Mq models are clearly intermediate cases, but it 
is not clear at this point if there is a distinct sequence or variation around a trend. Comparison 
of the Ye maps to the velocity maps show that the vast majority of material deviating from the 
Ye ^0.5 background is moving outward, and for the case of the 12 Mq model remains so after 
an additional 400 ms of evolution. For all models we find both lower- and higher- ejecta with 
complex, intermixed distributions. Because detailed nucleosynthesis must await the completion of 
the explosive phase before a definitive measure of the composition of the ejecta can be made using 
our tracers, we refrain from any premature quantification of the ejecta Ye except to highlight the 
potential for considerable differences in ejecta composition across progenitor mass. 

4. Summary and a look ahead 

Our four axisymmetric models all show large scale shock expansion to many thousand km with 
'diagnostic' energies of significant fraction of the canonical 1 Bethe(= 10^^ erg) explosion energy 
[5]. Though the 12 Mq model at 850 ms may have completed the accretion/neutrino heating phase, 
the other models at only 500 ms are clearly incomplete. From the models we can view a number 
of primary and secondary quantities to better understand how these explosions occur. The pattern 
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Electron fraction; 15 Solar mass model; 450 ms 
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Figure 3: Electron fraction, 7^, at 450 ms after bounce for all models. 

of neutrino heating and cooling can show the source of initial neutrino-driven convection and the 
intense heating at the bottom of the accretion funnels formed by the curvature of the aspherical 
shock. The flow of the accreted fluid into these funnels can also be clearly seen by following our 
tracer particles, which will eventually be used to compute detailed nucleosynthesis of the ejecta in 
our finalized simulations. 

The imposition of axisymmetry can affect the development of CCSNe and to understand these 
effects requires models without imposed dimensional symmetries. We are working to integrate 
the 'Yin- Yang' grid [19] into Chimera. The 'Yin- Yang' grid gains advantage by eliminating 
the constricted polar zones of standard spherical-polar grids and the time-step limitations that they 
bring. During our testing phase, we constructed a low-resolution (^33°) pilot run, initiated from 
a perturbed spherically symmetric model ^50 ms after bounce. This pilot run, over the course 
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of 100 ms evolution, shows neutrino driven convection as well as hints of an eariy spiral SASI 
mode. The shock position (viewed in a slice through the center) appears to rotate Uke a disk 
with an off-center pivot on to which is superimposed the distortions of of the shock surface from 
rising convective plumes. We have recently initiated a better resolved (^1.3°) run using the same 
complete physics set of our 2D simulations. We will present details of this model in future papers. 
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